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Abstract 

It is known experimentally that an imperfectly 
expanded rectangular jet from a thin-lip convergent 
nozzle emits only a single dominant screech tone. 
The frequency of the screech tone decreases contin- 
uously with increase in jet Mach number. However, 
for a supersonic jet issued from a bevelled nozzle 
or a convergent-divergent nozzle with straight side 
walls, the shock cell structure and the screech fre- 
quency pattern are fairly complicated and have not 
been predicted before. In this paper, it is shown that 
the shock cell structures of these jets can be decom- 
posed into waveguide modes of the jet flow. The 
screech frequencies are related to the higher-order 
waveguide modes following the weakest-link screech 
tone theory. The measured screech frequencies are 
found to compare well with the predicted screech 
frequency curves. 

1. Introduction 

Recently, there is a renewed interest in the pos- 
sible application of rectangular jets to future gen- 
eration aircraft propulsive system. Rectangular jets 
are regarded to have the potential for jet mixing en- 
hancement and noise reduction. Furthermore, they 
are known to have thrust vectoring capability. 

When operating at imperfectly expanded con- 
ditions, a supersonic rectangular jet emits discrete 
frequency screech tones. Earlier, these tones have 
been studied experimentally by Powell 1 , Hammitt 2 , 
Krothapalli et a/. 3 and others. During the last 
few years, they are the subject of several research 
studies 4 " 10 . In the most recent work of Raman 10 , 
bevelled rectangular nozzles were used. By shaping 
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the nozzle exit geometry, it was hoped that there 
could be additional noise and mixing benefits. 

The shock cell structure of an imperfectly ex- 
panded rectangular jet with a straight exit is quite 
complicated. Figure 1 shows the spark schlieren im- 
age of the top view of such a shock cell structure 
associated with an aspect ratio 5 convergent rect- 
angular nozzle. The jet Mach number is 1.4. If, 
instead of a regular straight exit nozzle, a bevelled 
nozzle is used, the shock cell structure becomes even 
more complex. This can be seen in figure 2. The jet 
Mach number is again 1.4. The shock cell structure 
is not only highly three-dimensional but also con- 
tains many fine-scale features. However, overall the 
shock cell structure is still quasi-periodic. 

It was observed experimentally, since the pio- 
neering work of Powell 1 , that a rectangular jet from 
a thin-lip convergent nozzle emits only one dominant 
screech tone (and its harmonics). Over the years, 
this finding has been confirmed by all subsequent in- 
vestigators. Figure 3 shows a typical variation of the 
screech frequency with increase in jet Mach number. 
The data 10 is for a jet issued from an aspect ratio 5 
convergent nozzle with exit dimensions 6.58 cm by 
1.32 cm. In contrast to the simplicity of the reg- 
ular straight nozzle case, the tone frequency Mach 
number relation for rectangular jets from bevelled 
nozzles, is highly complicated 10 . Observations re- 
vealed that it is possible for a jet to emit more than 
one dominant screech tones simultaneously. A closer 
inspection of the data indicates that even sudden fre- 
quency jumps can occur. At the present time, the 
physical mechanism responsible for this screech tone 
frequency pattern is not well understood. 

The primary objective of the present investiga- 
tion is to seek a physical explanation of the observed 
screech frequency jet Mach number relation for bev- 
elled rectangular nozzles. Both convergent and 
convergent-divergent nozzles are considered. Obvi- 
ously, the screech frequency must be closely related 
to the shock cell structure in the jet plume. For this 
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reason, a linear model shock cell structure of a super- 
sonic jet from a bevelled rectangular nozzle will first 
be established. The use of a linear shock cell model 
for screech frequency and broadband shock associ- 
ated noise prediction is not new. Such a model has 
been used successfully in ref. [11] to [15]. As in pre- 
vious works, the present shock cell structure jnodel 
is formed by a superposition of the waveguide modes 
of the jet flow. For screech frequency prediction, at- 
tention is focussed on the lower-order modes. 

It is known 1,16 that screech tones are generated 
by a feedback loop. According to the weakest-link 
theory 17 , the feedback loop is driven by the large- 
scale instability waves of the jet flow. The instability 
waves are excited at the nozzle exit region near the 
nozzle lip, where the jet is most receptive to external 
excitation, by acoustic disturbances. The excited in- 
stability waves extract energy from the mean flow 
and grow quickly as they propagate downstream. 
After propagating a distance of about three to four 
shock cells, they, having grown to relatively high am- 
plitudes, interact strongly with the shock cell struc- 
ture. As a result of the interaction, strong acous- 
tic waves are generated. Part of the acoustic waves 
propagate upstream outside the jet. Upon imping- 
ing on the jet near the nozzle lip region, the acoustic 
disturbances excite the shear layer, thus generating 
new instability waves. In this way, the feedback loop 
is closed. Tam, Seiner and Yu 17 observed that the 
acoustic waves that radiated to the nozzle lip re- 
gion were confined to a narrow frequency band if 
they were generated by the interaction between the 
large-scale instability waves and the quasi-periodic 
shock cell structure of a supersonic jet. The band is 
centered around frequency f p given by, 

f - ^ (1) 

h 2*1 + £) U 

where u c is the convection or phase velocity of the 
instability waves, is the ambient sound speed 

and k is the wavenumber of the shock cell structure. 
The weakest link of the feedback loop is the pro- 
cess of excitation of the large-scale instability waves 
by the feedback acoustic waves near the nozzle lip. 
In order to maintain the loop, the screech tone fre- 
quency must fall within the narrow band centered 
around f p given by (1). Because of this, equation 
(1) offers a good estimate of the screech frequency. 

To be able to use (1) for screech frequency pre- 
diction, u c and k must first be found. For cir- 
cular jets, Harper-Bourne and Fisher 18 had mea- 
sured experimentally that u e — O.luj ( Uj is the 
jet speed). However, for rectangular jets Tam and 


Reddy 19 found it was more appropriate to use a 
lower convection velocity. They recommended 


u c ~ O.bhuj. 


( 2 ) 


In this work, formula (2) will be used throughout. 

By means of a vortex sheet jet model, Tam 1 ” 
found that the pressure distribution, p, associated 
with the shock cell structure of a rectangular jet can 
be decomposed into waveguide modes with eigen- 
functions (mode shape), <i>nm{y, z ) and wavenumber 
k nm . With respect to a Cartesian coordinate sys- 
tem shown in figure 4 for a jet with fully expanded 
width bj and height hj , p may be represented by an 
expansion in the form 

oo cc 

p( x i V) z ) — } ] ^nmiy 1 *)[Anm COS (k nm x) 

n = l m— 1 

+ Bnm sin {k nm x)] 


where 



(4) 

(5) 


The coefficients A nm and B nm are to be determined 
by the nozzle exit conditions, bj and hj are related 
to the nozzle exit dimensions, 6 and h, by 


where 
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h 
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\bh )\b + h , 
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(7) 
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> + l 

(8) 




In the above, Mj is the fully expanded jet Mach 
number and Md is the nozzle design Mach number. 

For shock cell mode prediction, the linear model 
above is fairly accurate. For instance, the screech 
frequency of the measured data in figure 3 can be 
predicted by taking k in formula (I) to be the 
wavenumber of the lowest waveguide mode (n = 1, 
m = 1); i.e., 
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the initial conditions is 


The calculated values form the solid curve of this 
figure. As can be seen, there is good agreement be- 
tween the measured data and the theoretical predic- 
tion over a wide range of Mach numbers. 

For bevelled rectangular nozzles, the associated 
shock cells are necessarily more complicated. For 
the stated purpose of this investigation, namely, to 
obtain a physical understanding of the screech fre- 
quency pattern of these jets, it is believed that the 
use of a linear shock cell model is sufficient. The de- 
velopment of such a model will be discussed below. 
The results of the model is then used in section 3 
of this paper to predict the screech tone frequency 
jet Mach number relation. Comparisons with ex- 
perimental measurements are carried out. The ex- 
perimental measurements involved jets from single 
bevelled nozzles, double bevelled nozzles and regular 
nozzles. Both convergent and convergent-divergent 
nozzles are included. 

2. The Shock Cell Structure and 
Screech Frequencies 

2.1 The Physical Model 

Let us consider the shock cell structure of a rect- 
angular jet from a single bevelled rectangular nozzle. 
Figure 5 shows the plan view of such a jet. For sim- 
plicity, we will use a vortex sheet jet model as in the 
work of Tam 12 . We will assume that the jet Mach 
number Mj is such that the bevelled angle /? is larger 
than the Mach angle p. When this is true the distur- 
bances generated by the upstream nozzle side wall 
would not reach the inside of the nozzle. That is, 
all the flow disturbances produced by the boundary 
discontinuities at the nozzle exit are confined within 
the jet flow. 

For weak shock cells, the shock structure dis- 
turbances are governed by the linearized equations 
of motion. On following Tam 12 , it is straightforward 
to find that the pressure, p, of the shock cells satisfies 
the equation 


V 2 P-M/|^ = °. (10) 

On the boundaries of the jet (see figure 5) the pres- 
sure balance boundary condition across the vortex 
sheet must be satisfied. That is, at y — 0, y = tj, 
z = 0, z = hj 

P = 0. (11) 

At the nozzle exit, y = — x tan/?, both p and are 
prescribed. These are the initial conditions deter- 
mined by the supersonic flow upstream (inside the 
nozzle). For a convergent nozzle, a simple model of 


p = Ap 



(constant) 


(12) 

(13) 


For convergent-divergent nozzles with straight side 
walls, Tam and Reddy 19 have pointed out that a 
second set of shock cells would be generated at the 
nozzle throat. Thus the pressure distribution at the 
nozzle exit cannot be constant. Therefore, (12) and 
(13) would not be appropriate for these types of noz- 
zles. 

2.2 Numerical Solution 

Initial boundary value problem (10) to (13) can 
be solved relatively easily if the solution is first ex- 
panded as a Fourier series in z. That is, let 

P(x,y,z) = p m (i,y)sin • (14) 

On substitution of (14) into (10), the governing 
equation for p m (x,y) is found to be, 




5 0 ^ 0 0 

Pm , m ~ 7r ~ ~ 
dy 2 h) Pm 


= 0 . 


(15) 


The boundary conditions from (11) are, at y = 0, 
V = bj 

Pm = 0. (16) 

On expanding (12) and (13) in Fourier series of 
sin(“ L£ ), it is straightforward to find that the initial 
conditions at the nozzle exit, y = — x tan /?, are 

2 

p m = (1 - cos nu)Ap 

rmr 

^ = o. 

dx 

We solve equation (15), boundary condition 
(16) and initial conditions (17) and (18) by finite 
difference method. To keep the order of the finite 
difference equations the same as the order of the 
original partial differential equation (15), we choose 
to approximate the spatial derivatives by second- 
order central differences. One principal reason for 
using a second order finite difference approximation 
instead of a higher-order scheme is that the number 
of initial and boundary conditions for the resulting 
finite difference equation is identical to those of the 
original partial differential equation problem. Had 
a higher-order scheme been used, extra initial and 
boundary conditions must be imposed in order to 
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obtain a unique solution. Then a way must be found 
to impose these extra conditions so that they have 
no effect on the numerical results. With respect to 
a mesh of spacing Ax and Ay as shown in figure 6, 
the discretized form of (15) is, 


, ,^2 IN (Pz+l.fc 

( m j ~ i ) ^ 

_ (Pf,fc + 1 ~ 2p<,fc +_Pt_ 1 k- 1) 
A y 2 


(19) 


m 2 7r 2 

r 


Pt,k = ° 


where (£, k) are the mesh indices in the x- and y- 
directions. (19) may be rewritten in the following 
form so that the value of p t + l k may be computed 
once its values at the last two columns, £ and £ — 1, 
are known. 



Suppose Ay = ^ where K is a large integer. Then 
(20) is valid for 1 < k < K — 1. 

Boundary conditions (16) yield, 

Pt , o — Pi,k ~ 0- (21) 

To ensure that the computed solution satisfies initial 
conditions (17) and (18) at the nozzle exit, a simple 
way is to start calculating the solution using (20) at 
the column AA ' {£ = —N) as shown in figure 6. We 
set 

2 

p_w.it = P-(jv+i),fc = — (1 - cos m7r)Ap (22) 

and use (20) to march the solution in the x-direction. 
As soon as the value of p t k at a new £ column (at 
the first step, this will be £ = — N + 1) is found all 
the values of k at mesh points that lie inside the 
nozzle are reset to ~:(1 ~ cost )Ap. This assures 
that at the nozzle exit the slope in the x-direction 
is zero and also the pressure there is equal to the 
right side of (17). The marching process is to be 
continued until £ > 0 is reached. 

Finite difference equation (20) is computation- 
ally unstable unless a small value of Ax is used. The 


stability requirement of (20) can be established by 
performing a standard Von Neumann analysis. Let 

P<, t = *V At . (23) 

Substitution of (23) into (20) leads to the following 
equation for a. 


cr 2 — Qa +1—0 (24) 


where 






(cos A — 1) 


m 2 T 2 (Ax) 2 

+ h) Mj - 1 


The roots of (24) are 

<r± = \[Q±(Q :i -^}. 


(25) 


(26) 


One of the roots of (26) will have a magnitude 
greater than unity unless <r± are complex conjugates 
of each other in which case |<r±| = 1. The compu- 
tation scheme is, therefore, stable if Q~ < 4. This 
yields 


-4 < 




( cos A - 1 ) 


m 2 T 2 (Ax) 2 
k] Mj - 1 


(27) 


The right side inequality of (27) is always satisfied. 
The left side inequality would be satisfied for all A 
provided, 


/Ax\ 2 4 m 2 T 2 (Ax) 2 

\ Ay/ + h) Mj - 1 


< 4. 


Thus, numerical stability is assured if Ax is taken 
to be 

(Mj - 1) = Ay 


Ax < 


[l + m^i(M)2]i 


(28) 


In the present numerical computation, the step size 
Ax used is always smaller than that given by the 
right side of (28). 


2.3 Waveguide Modes of the Shock Cell Structure 

In the region of the jet downstream of the y- 
axis at £ = 0, the shock cell solution can again be 
expanded in terms of the waveguide modes of the jet 
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flow as in equation (3). In this case p m (x , y) may be 
written in the form, 


CO / \ 

p m { x >y) = X! sin 


n = l 


(29) 


+ B nm s\n(k nm x)]. 


The coefficients A nm and B nm are to be determined 
by joining expansion (29) smoothly to the numeri- 
cal solution at a location, say, at £ = L (L > 0). 
The joining conditions are the continuity of p m and 

On replacing integrals by appropriate sums 
(the numerical solution is in discretized form), it is 
straightforward to find from the joining conditions 
that the values of A nm and B nm are given by, 


Anm — 


2A y 

knm^j 


K 


knm COs(t nm LAx) ^ ^ PL,k 
rnrkAy N 


k-Q 

Sin(^r»m-^ / Ax) 


K 


^ 9A X (?£, + !,*: ~ PL-\,k) 


k=0 


( rnrkAy \ 


(30) 


Bnm 


2Ay 

krim^j 


K 


knm sin(fc nm LAx) ^ ^ PL,k 


kzzO 


( rnrkAy\ 


• sin 


K 




+ cos(fc„ m LAx) 


1 

^ 2Ax^ L + 1,Jk ~PL-\,k) 


<fe = 0 


( rink Ay \ 

~JJ 


(31) 


The amplitude of the (m, n) waveguide mode is 


C nm = {A 2 nm + B 3 nm )i. (32) 


2.4 Screech Frequencies 

Beyond the most downstream edge of the noz- 
zle, the shock cell structure is quasi-periodic. This is 
the region where the feedback acoustic waves of the 
screech feedback loop are generated. The physical 
situation is very similar to that of a jet with a regu- 
lar straight nozzle exit. The major difference lies in 
that the dominant waveguide modes that make up 


the quasi-periodic shock cell structure are not nec- 
essarily the lowest order modes. Moreover, it is pos- 
sible that more than one mode is dominant. As far 
as the screech tone feedback loop is concerned, the 
weakest-link theory is still applicable. However, we 
expect the resulting screech tones to be more com- 
plex. This is because they are related to the higher- 
order modes as well as more than one feedback loop 
may be operating at the same time. On following 
Tam 12 , the screech frequencies can be calculated by 
replacing k in equation (1) by k nm given by equa- 
tion (5). For a large aspect ratio jet (6 ; >> hj) only 
the m = 1 modes are relevant. On restricting our 
consideration to the m = 1 modes, we obtain, 


/n = 



(.'// - n* 


(33) 


n = 1,2,3,.... 

3. Comparisons with Experiments 

Detailed variations of the screech tone frequen- 
cies with jet Mach number for an aspect ratio 5 jet 
have been measured by Raman 10 . In this section, we 
will compare the frequencies of his measured data 
with formula (33). 

In the experiments of Raman, six rectangular 
nozzles were used. Three were convergent nozzles 
and the other three were C-D nozzles. In each fam- 
ily of nozzles, one has a regular straight exit ge- 
ometry, the second has a single-bevelled and the 
third a double-bevelled geometry. The dimensions 
of these nozzles were given in the works of Rice and 
Raman 5,6 . The variation of the screech tone frequen- 
cies with fully expanded jet Mach number for each 
nozzle was measured. However, because more effort 
was spent on the C-D nozzles, only data related to 
them were reported. In addition to screech frequen- 
cies, the spanwise phase difference and transverse 
phase difference of the screech tones were also mea- 
sured. According to the phase information, Raman 
was able to separate the screech tones into different 
modes. 

Figure 7 shows the measured screech tone fre- 
quencies as a function of the jet Mach number tor 
the supersonic jet from the single-bevelled conver- 
gent rectangular nozzle. The data does not seem to 
form a continuous curve. There are frequency jumps 
at several jet Mach number. Shown also in this fig- 
ure are the screech frequencies of formula (33) with 
spanwise mode number n = 1,2,3, and 4. As can 
be seen, the data correlate well with the theoretical 
frequency curves over the entire range of measured 
Mach numbers. Figure 8 shows a comparison of the 
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measured screech frequency and the frequencies of 
formula (33) for rectangular jets from the double- 
bevelled convergent nozzle. The screech frequencies 
agree well with the n = 5 and n = 6 frequency 
curves. 

Figure 9 shows the variation of the screech fre- 
quency with Mach number for the supersonic jet 
from the C-D rectangular nozzle with a regular exit 
geometry. Based on his spanwise phase data, Raman 
was able to identify three distinct screech modes. 
Their frequencies fit very well with the frequencies 
of the n = 1,3 and 5 modes of formula (33). Fig- 
ure 10 provides the measured screech frequency jet 
Mach number relation of the supersonic jet from the 
single-bevelled C-D rectangular nozzle. According 
to the phase data, there are three modes, i.e., mode 
I, modes III A and IIIB. Shown also are the screech 
frequency curves for n = 1,2,3 and 4 modes ac- 
cording to formula (33). The data appear to agree 
well with the frequency curves of the n = 1,2 and 4 
modes. Finally, figure 8 shows the measured screech 
frequency variation with jet Mach number of the 
supersonic jet from the C-D double-bevelled rect- 
angular nozzle. Two modes have been identified. 
Their frequencies are in reasonable agreement with 
the n = 6 and n = 7 frequency curves of formula 
(33). 

Formula (33) predicts many screech frequency 
bands. However, not all the tones of the frequency 
bands are excited. We believe that the selection of 
the excited frequencies is determined largely by two 
factors. They are the spatial growth rates of the in- 
stability waves of the jet flow and the strength of the 
waveguide modes of the shock cell structure. Here, 
we will examine, by means of the shock cell model 
above, how the strength of the waveguide modes in- 
fluences tone selection as well as the frequency jumps 
in the data of figure 7. For this purpose, the am- 
plitudes of the waveguide modes corresponding to 
m = 1, n = 1,2, 3, 4 were calculated by the march- 
ing scheme (20). In the numerical computation K 
was taken to be 8000. The size of the marching step 
Ax was chosen to be smaller than that given by (28). 
The solution was marched from the starting line AA! 
in figure 6 to a line downstream of the y-axis before 
the coefficients A„i, B n \ and C n \ were evaluated ac- 
cording to equations (30) to (32). To check the accu- 
racy of the numerical results, each computation was 
marched further downstream. At this new location, 
the amplitude C n \ was again determined. It was 
found in every case that the two calculated values 
of Cni agreed to, at least, three significant figures. 
This offers assurance that the numerical results are 


accurate. 

Figure 12 shows the computed waveguide mode 
amplitude C n i, n = 1,2,3, 4 of the shock cell struc- 
ture for a single bevelled convergent rectangular noz- 
zle used in the experiment of Raman 10 . In the com- 
putation, A p was set equal to unity for convenience. 
At low supersonic Mach number, the dominant mode 
is the n = 4 mode. As the Mach number increases, 
the dominant mode shifts to n = 3. then n = 2, 
and finally at Mj > 2 the ti = 1 mode becomes the 
most dominant. This trend correlates well with the 
measured data in figure 7. 

If the shock cell strength is the only factor that 
determines the dominant screech tone, then the in- 
formation in figure 8 can be used to predict the lo- 
cations of the frequency jumps observed in figure 7. 
Let us assume that a shift in the excited frequency 
is possible when the waveguide mode amplitude of 
the adjacent band is within 15% of each other. Fig- 
ure 13 shows the mode switching Mach number or 
the expected locations of frequency jumps. The pre- 
dicted locations for mode switching from n = 4 to 
n = 3 and from n = 3 to n = 2 (indicated by ar- 
rows) are in fair agreement with measurements. The 
switch from n = 2 to n = 1 mode takes place around 
Mj = 1.75 whereas the calculated value falls around 
Mj = 2.0. We believe that the discrepancy is solely 
because the effect of the growth rate of the insta- 
bility waves has not been taken into account. De- 
spite this discrepancy, the good overall agreement 
leaves little doubt that the strength of the waveguide 
modes of the shock cell structure plays a crucial role 
in the selection of the observed screech tones. 

4. Summary and Discussion 

The screech tone frequency patterns of super- 
sonic jets from bevelled rectangular nozzles and C- 
D nozzles with straight side walls are quite complex. 
In this paper, it is suggested that these tones are 
related to the higher order waveguide modes of the 
shock cell structure in the jet plume. By using a 
simple vortex sheet jet model to calculate the linear 
shock cell waveguide modes of the jet and then calcu- 
late the screech frequencies by the axial wavenumber 
of these modes, it is shown that the predicted screech 
frequency curves agree well with all the measured 
data. The good agreement provides strong support 
for the waveguide mode concept of the shock cell 
structure and the weakest-link feedback theory of 
jet screech. 

In the present work, the shock cell structure 
is modelled by a superposition of waveguide modes 
of the jet flow. This modal decomposition is un- 
questionably at odds with the concept of shock cells 
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with sharp discontinuities. One major difference be- 
tween the two models is that there can only be one 
screech tone associated with a shock cell structure 
with sharp discontinuities (there is only one shock 
cell length). On the other hand, there can be more 
than one screech tone according to the waveguide 
model. In the above, we have found good ^agree- 
ments between the screech frequencies predicted by 
the waveguide shock cell model and experimental 
measurements. Furthermore, Raman 10 has shown 
experimentally that two screech tones can exist si- 
multaneously. Thus there is no question that a 
waveguide model representation of the shock cells 
of a supersonic jet is valid. However, the authors 
believe that the range of validity is limited only to 
the moderately imperfectly expanded Mach number 
range of the nozzle. For highly imperfectly expanded 
jets, nonlinear effects dominate the formation of the 
shock cell structure. In this case, the structure is 
made up of sharp shocks. It follows, therefore, for 
highly imperfectly expanded jets there would only 
be a single screech tone. This appears to be consis- 
tent with all the screech tone data available in the 
literature at the present time. 
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Figure 1 Spark schlieren image of the shock cell structure 
of a Mach 1.4, aspect ratio 5, rectangular jet from a regular 
convergent nozzle 



Figure 2 Spark schlieren image of the shock cell structure 
of a Mach 1.4, aspect ratio 5, rectangular jet from a single 
bevelled convergent nozzle 




Figure 3 Variation of the screech frequency with fully expanded 
jet Mach number for a supersonic jet from a convergent rectan- 
gular nozzle. Aspect ratio = 5. Exit geometry = regular. O 
measurements 10 ; theory 1 *. 
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Figure 6 The computation domain with mesh size Ax and Ay . 
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Figure 4 Coordinates used to calculate the shock cell structure of 
rectangular jets. 



Figure 7 Variation of the screech frequencies with fully expanded 
jet Mach number for a supersonic jet from a convergent rectangular 
nozzle. Aspect ratio = 5. Exit geometry = single bevelled. A 

measurements 10 . Equation (33); n=l. n=2, — • — 

n=3, n=4. 
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Figure 5 Top view of a vortex sheet jet from a single bevelled 
rectangular nozzle. 


Figure 8 Variation of the screech frequencies with fully expanded 
jet Mach number for a supersonic jet from a convergent rectangular 
nozzle. Aspect ratio = 5. Exit geometry = double bevelled. □ 

measurements 10 . Equation (33); n = l, n— 2, — - — n=3, 

n=4, — •• — n=5, n=6. 
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Figure 9 Variation of the screech frequencies with fully expanded 
jet Mach number for a supersonic jet from a convergent-divergent 
rectangular nozzle. Aspect ratio = 5. A/j = 1.4. Exit geometry 
= regular. Measurements 10 : 0 mode 1-1, A mode 1-2, □ mode II 
. Equation(33); n=l. n=3 t — • — n=5. 


Figure 11 Variation of the screech frequencies with fully expanded 
jet Mach number for a supersonic jet from a convergent-divergent 
rectangular nozzle. Aspect ratio = 5. Exit geometry = double 
bevelled. Measurements 10 : O mode I. A mode II. Equation(33): 
n=6, n=7. 




Figure 10 Variation of the screech frequencies with fully expanded Figure 12 The amplitude of the waveguide modes of the shock 
jet Mach number for a supersonic jet from a convergent-divergent cell structure in the plume of a supersonic jet downstream of a 
rectangular nozzle. Aspect ratio = 5. = 1.4. Exit geometry convergent single bevelled rectangular nozzle. Aspect ratio = 5. 

= single bevelled. Measurements 10 : O mode I, A mode III- A. □ Cii * ' * ^ 31 • 

mode III-B . Equation(33}; n=l, n=2, — • — n=3, — 

n=4. 







Figure 13 Mode-switching Mach numbers, llll based on com- 
puted waveguide mode amplitudes. 



